GE 


^vMM  ■>««  — 

»«< 

taHinwii  ai  m^Dmxw-. . 

]im-4N)  w4«oi«omc»«i  <»— 


1.  A«iNOr  USi  OW.Y  a««v«  (H4nii)  1 1.  ■! HMT  OATI 


4.  mu  ANO  sutmu 

On  The  stability  of  Finite  Elastic  Deformations 


om  m.  0A»A-or«t 


W">*.  ow  «Mt  t«f  itmr^merm.  itwtXM  •(•km  mm  to^tm, 

mutaaimm 

mtanan  wwtwL  Otunartw  hr  mt^rmanm  QM»»iMrii>  mt  um^n. 
tliwriiow  %#>nioii  Pr^taa  (8;»*Ai«ll  (X  Wl 


gan  Technological  Univ. 
.on,  MI  49931 


I.  SMMSOlUMe/MONITOAiNa  ACCNOr  NAMI(S)  ANO 

U.  S.  Araj  iUscarch  Office 
P.  0.  Box  12211 

Reacarcli  TrUngle  Park,  NC  27709-2211 


S.  rUN0<NQ  NUMCfAS 


t.  OKCANIZATION 

Rci^oaT  NUMia 


10.  SRONSOUNC/  MONiTOMiNO 
ACfNCr  UKMT  NUMtCR 


/iiU)  "^6- 


11.  SUmiMCNTARV  NOTIS 

The  vlav,  oplolona  and/or  flndinga  contalnad  In  this  report  are  those  of  the 
author  (a)  and  should  not  be  construed  as  an  official  Department  of  the  Amy 
position,  policy,  or  decision,  unless  so  designated  by  other  docuauuntatlon. 


12«.  OnTRISUTtON/AVAIlASajTT  STATIMINT  12b.  OtSTIUSUTKM  COOC 

Approved  for  public  release;  dletrlbutiom  unlialted. 


1A  AS5TRACT  (kUidmumiOOwofta^ 

The  stability  of  inelastic  deformation  and  the  formation  of  shear  bands  are 
examined  by  considering  a  recently  developed  gradient-dependent  theory  of  elasto- 
viscoplasticity.  The  main  is  ues  addressed  in  the  present  work  are  the  role  of 
viscosity,  inertia  and  higher  order  strain  gradients  in  providing  internal  length 
scales,  as  well  as  the  effect  of  back  stress  (deformation  induced  anistropy)  and 
plastic  spin  (texture  development)  on  the  structure  of  the  shear  bands  in  the  post 
localization  regime. 


93-03326 


llllllllii! 


14.  SUUta  TtWMS 

stability  of  Plastic  Deformations 


IS.  NUMtfR  Of  AAGfS 


14.  naa  cooc 


1>.  SICUWTV  OASStflCATION  I  It.  SICUWTY  CLASSIfICATJON  I U.  StCUIIfTY  OASStftCATlOM  i  20.  UMITATION  Of  AISTRAO 

Of  RIKMT  I  or  THIS  PAGI  I  Of  AtSTlua  I 


OKCLASSIFIED 


MSN  7540-01 -280.5500 


UNCLASSIFIED 


UNCLASSmSO 


StAodard  form  299  (Aav  2  89) 
Pmarn4  I,  M«i  Ma  m-'i 
m-toa 


On  the  Stability  of  Finite  Plastic  Deformations 


H.M.  ZWb 

Washington  State  Umversty 
Pullman,  WA  99164-2920.  USA.  and 

E.C.  Aifantls 

Michigan  Technotogiot  University 
Hougton,  Ml  49931 .  USA  and 
ArtstoSe  University  ol  Thessaloniki 
Thessakxilki  54006,  GREECE 


ACCc  :0  1  To' 


MiS 

D  i  i ' .  \  j-u 

U, .  r:  T-  i  i 


7 


By 

Diit  ! 


Av.:riub;!ity  Co.UjS  \ 

.....  1 

Dist 

Al 

Avail  . 

:  Of 

Cldl 

Summary 

The  stability  at  ineiastc  detormatlon  and  the  lormaiion  oi  shear  bands  are  enamined  by  consldedno  a 
recently  developed  gratfent-dependent  theory  ol  elasto-viscopiastidty.  The  main  Issues  addressed  In 
the  present  work  are  the  role  ot  viscosity.  Inertia  and  higher  order  strain  gradtonts  In  providing  internal 
len^h  scales,  as  we«  as  the  etfect  of  back  stress  (deformation  Induced  anisotropy)  and  plastic  spin 
(texture  development)  on  the  structure  of  the  sheer  bands  in  the  post  locsizatlon  regime. 


Introduction 

The  problems  ol  consinutive  modelling  ol  finite  plastic  deformations  and  the  associated  ohenomera  of 
deformation  pattom-tormafioo  and  shear  banding  have  been  the  subjects  ol  Interest  among  m»iy 
investigators  during  trie  past  two  decades;  recent  reviews  can  be  found  m  the  works  of  Rice  J 1  ],  Ailwas 
12a.b)  and  Zbib  and  Aitartis  (Sa.b).  When  examining  these  phenomena,  theoudai  questions  that  arise 
are  the  identilication  oi  the  basic  scaling  mechanisms,  the  associated  microsbuctural  aspects  (see  lor 
example,  Korbel  et  al.  (4().  the  appropriate  internal  variables  that  a  oonsttiuti^  theory  should  inckato 
and  the  way  the  macroscopic  constitutive  response  Is  reined  to  the  physicai  microprocesses  (see  lor 
example.  AHantis  (2]  and  Orucksr  {5]). 

The  stability  of  plastic  low  and  the  shear  banding  phenomenon.  In  particulaf.  have  gained  the  attertion 
of  many  researchers  in  both  the  metallurgical  and  mechanics  communiiss  (1-1 1],  While  the  nain 
mathematicat  and  physical  aspects  pertaining  to  the  understanding  ol  the  mechanisms  leading  to  Vie 
onset  of  the  phenomerxin  have  been  examined,  work  on  the  charactenzaion  and  evolution  of  shear 
bands  including  the  determination  ol  band  widths  and  spadngs  has  only  recently  begun.  The  miin 
difliculty  in  this  problem  lies  m  the  proper  modelling  ot  material  behavior  in  tie  post-localization  reg^ 
once  the  shear  bands  are  formed  This  can  be  attibuted  to  the  (act  that  in  tvs  regime  the  deformation 
field  becomes  highly  iihomogenAous  and.  therefore,  constitutive  equations  originally  proposed  for 
homogeneous  or  ‘near  homogeneous’  fields  become  insufficient.  To  overcome  such  difficulties,  it  aw 
suggested  in  the  work  of  Aifantis  and  co-workers  p.3, 12-16]  to  include  the  higher  order  strain  gradieMs 
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Into  ttw  oonstKutive  theory  ol  plasticity  to  aooount  lor  the  heterogeneity  al  plastic  flow  and  capture 
Internal  detormation  wavelengths 

In  the  present  worir  we  first  review  the  'scafa-invarlanca*  approach  outSned  in  (21  arid  used  in  (3}  to 
obtain  physically  based  oorstitutive  aquatons  (or  the  plastic  strstcmng.  bade  sbess  (modelling 
anisotropy)  and  plastic  spin  (modelling  texiwe  development)  Such  a  theory  may  be  viewed  as  a 
compromise  between  *muitislip'  crystal  piasScily  averaging  models  and  pmty  continuum  plasticity 
models.  The  theory  s  then  modified  to  accotri  lor  the  heterogeneily  ol  piasbc  How  occurring  during 
the  shew  banding  process.  Thrs,  in  iact.  can  be  directly  related  to  iracrostnxlural  inhomoganeities 
resuldrig,  for  er^ample.  from  tangled  cell  structures  of  dislocations,  testuie  development,  etc.  This,  in 
turn,  leads  to  the  deveiopment  ol  strain  heteiogeneities  causing  the  intiatton  of  localUed  detonrabon 
bands.  The  macroscopic  manifestation  of  such  microscopic  everts  is  sssumed  here  to  be  the 
development  of  higher  order  strain  graOmta  which  would  enter  toto  the  form  of  macroecoplc 
phenomenological  theory.  In  fact  we  retain  tia  structure  of  the  dassicaf  dieory  of  plastidly.  Inebiding 
the  yield  surface  (lor  rate-independent  materMs).  the  normality  condWon  and  the  now  rule,  but  kibeduoe 
higher  order  strair>-gradieats  into  the  equatian  for  the  flow  stress. 

In  order  to  examine  toe  role  ol  higher  order  strain  gradients  on  toe  material  behavtor,  we  consider  a 
simpfified  orte-dimensiortal  probiem  arid  investigate  the  effect  ol  inerfla.  vtsoosily  and  gradients  on  the 
structure  of  (he  shear  band.  K  is  shown  toat  lor  the  quasi-static  loading  case,  higher  order  strain 
gradients  are  the  only  source  to  provide  m  internal  length  scale  to  (he  problem  lor  boto  rate- 
independent  and  rate-deperiderTi  materials.  For  the  dynamic  loading  case,  however,  It  i  shewn  toat 
viscosity  and  inertia  together,  are  sutl'icieni  to  prc'lde  a  natural  length  scale  to  the  problem  eitooul 
resorting  to  higher  order  strain  grarSents.  butimen  viscosity  is  dropped,  higher  order  gradients  become 
again  necessary 

In  order  to  examirw  the  efieci  of  plastic  spin  and  anisotropy  on  the  devetoprnent  of  shear  bands  we 
consider  a  two-dimenslortal  probiem.  The  problem  Is  treated  numerically  usirvg  ttve  tlnlie  tiemem 
method.  It  is  shown  that  for  too  case  ol  shear  banding  from  a  ground-stale  of  pure  exiensfon  (no 
shearing),  the  plastic  spin  has  no  influence  on  the  onset  of  instability,  but  It  does  affect  the  develapment 
ot  severe  localizaiion.  In  tact,  an  increase  to  the  degree  of  plastic  spin  decreases  'ductility,*  in  the 
sense  (hat  the  value  of  average  strain  to  senore  locali2ation  and  'toad  collapse'  decreases. 

A  Gradient- Deoendert  Theory  of  Etasto-Plasflcity 

We  consider  a  continuous  medium  under  dyrwnic  loading  coridllions  whose  motion  Is  governed  by  the 
linear  momentum  equation. 


17 


^  a  «  p  w,  0) 

nytiere  o  Is  the  Cauctiy  strsss  tensor,  p  is  tf»  mass  rlemlty,  and  v  «  Nw  retodty  vector  fteW, 
Assuming  that  the  solid  undergoes  large  detormattons,  we  can  decompose  #»  delormation  gradient  f 
Into  elastic  F*.  plastic  F*  and  rotation  R“  parts  such  Itra*  F  »  R"  P  F*  as  proposed  by  Zblb  and  Aitanbs 
{3a).  Following  their  rnethod.  one  can  then  obtain  me  decomposition 

0  •  0'  ♦  0'  .  <*  •  W  -  W'. 

where  the  stretching  0  is  the  symmetric  part  of  the  velodty  gradient,  07  and  O'  are  the  elastic  and 
plastic  stretching  respectively.  o>  is  the  spin  of  me  maieilai  (substrucaire),  W  ts  the  antisymmetric  pan 
of  the  velocity  gradtant  (spin  of  the  continuum)  and  VF  Is  the  plastic  spin.  The  elastic  strain  rate  0*  is 
given  by  Hooke's  law.  modified  for  large  deformations,  as 

D*  ’  (Cl  ’i  ;  i  -  4-»8‘S«*, 

where  (C*I  is  the  etastldty  tensor  and  &  is  the  oorotatlor^ai  stress  rate. 

Using  a  scale-invartance  argument  and  a  meximizatlon  procedure,  AHanib  and  co-wortcars  (e.g.  2,3] 
have  been  able  to  derive  explicit  expressions  tor  me  piasttc  stretching  Xf  and  plastic  spin  WF  based  on 
the  process  of  crystal  slip  artd  dislocation  glide.  For  kinematic  hardening  plasticity  models  a  turns  out 
that  the  plastic  anisotropy  can  be  represented  by  the  back  stress  a  xmosa  evolution  Is.  In  taa. 
determined  by  (he  plastic  spin  WF.  As  a  special  case  o(  this  (ormuWton  or«  obtains  the  lolowing 
constitutive  eguations 


0^  *  ^(S'-«0.  4  *  /O'-cya;  a  »  a-ue*ae;  (4^ 

««  •  W-W',  W'  .  C(a0'-0  V 

where  prime  irxScalas  deviatortc  while.  f>,  c  and  C  are  material  parameters.  (As  usual.  t> 

(2  0' .  iF)'*  denoies  the  etfective  strain  and  T  » (S' .  SVZ)”*  the  effective  stress]. 

Equation  (4),  suggest  that  the  plastic  spin  is  a  measure  of  nonooaxlaMy  between  the  back  streas  and 
the  stretching.  Its  implication  to  large  detormtSon  has  been  examined  in  detal  by  the  authors  {3).  In 
.Mirtlcular.  It  was  shown  that  tsF  has  a  oonsiderabla  eltect  on  the  development  of  axial  stress  (or  strain) 
in  a  fixed  (or  free)  and  torsion  of  a  cytindrtcal  bar  and  mat  the  value  ot  the  parameter  i;  can  be  obtained 
from  such  torsion  tesla.  (The  results  of  Zblb  and  ASantls  (3]  suggested  that  ^  has  very  small  values; 
in  lact,  by  detining  5  ■  a/o,.  where  o,  is  the  initial  yield  stress,  they  obtained,  lor  example,  a  •  4.5  tor 
AL- 1 1 00).  The  implication  ot  the  plastic  spin  to  loc^uation  is  examined  in  the  last  section  ot  this  paper 
in  connection  with  me  numerical  analysis 
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ftate-Deoandent  Matertate 

Ttw  flow  stress  T  appearing  in  aquaton  (4),  and  its  dependence  on  strain  and  strain  rate  determines 
whether  the  materiat  is  rate-independent  or  rate-dependert.  In  addition  to  this  classical  interpretation 
ol  the  flow  stress.  It  was  proposed  h  (2.3)  that  higher  order  gradients  ot  micro  (disiocation  densities) 
or  macro  (strain)  variabtes  be  included  in  the  respective  oonstitutive  equation  in  order  to  capture  the 
details  ol  tMormalion  patterning  and  shear  banding  phenomena.  For  the  case  ol  macro  shear  bands, 
for  example,  the  following  gradiensdependent  expression  tor  the  flow  stress  was  proposed 

^  (5) 

^ t,  •  c,v*Y*c^  ivyi*, 

where  k  d.y  )  is  the  usual  homogsneous  stress  and  C,  and  C,  are  the  so-caned  gredtent 

coefficients.  This  simple  expression  provides  an  infemal  length  to  the  theory  ol  plasftcity  and,  hence, 
makes  possible  the  investigation  at  pattem-torming  instabWes  in  the  deformalton  tWd. 

Upon  combining  equations  (3)  and  (4)  with  (2)  we  obtain  the  following  gracflent-dependent  constitutive 
equation  lor  rate-deper«deni  elastic-viscoplastic  materials  with  Isotropicfldnemattc  hardening 

&  .  [C'P-^fC'KS'-.');  R  .  -  M, ; 

2i.  y  y 

where  4  can  be  defined  as  the  nonHnear  total  viscosity,  with  a  homogeneous  part  and  an 
inhomogeneous  part  pr 

When  elasttc  effects  are  neglected  9>*<0)  lor  targe  plastic  deformations,  ffien  we  obtain  Vie  following 
constitutive  relation  tor  a  gradient-dspendent  vlsooplastic  material  with  isobopicfldneniatic  hardening 

S  •  ZpOve-pt;  ^ 

where  p  «  -  tr3/3  is  the  hydrostafc  pressure.  This  equation  is  reminiscent  to  that  desalbing  a 
Newtonian  viscous  fluid  (without  a  and  a  gradient-dependent  viscosity). 

Rate- Independent  Materials 

The  dependence  of  the  flow  stress  r  on  J  indicafes  viscoplastic  behavior  Plastic  behavior  is  obtained 
by  dropping  itr  (5)  the  dependence  on  f  In  this  case  t  and  the  condition  of  loading  or  unloading  are 
obtained  from  the  yield  function  F  and  the  consistency  condition  F  -  0.  The  loading-unloading 
conditions  are  given  by 
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-  I  F  -0  and  r  *0  Of  1 
F’yJ,-T.  I  loading  or  nautrai  loating 
I  and  T  <0-  untoading 

whofo  4  »  {S'-a')-{S'  -«V2  oonsistoncy  condition  F-  0  Ifion  jdolds 
cy  f  .2C,7 ?  V  T  -  •  -<^''*1^^' ; 


wfiore  the  ‘altecilva*  hardenlr^g  modUus  h  is  givan  by 


e  '.(S'-aV^J^-^V^-^  1^1*. 

dt  2  2t  dy  dT  <*t 


(8) 


(9) 


(10) 


Relation  (9)  is  a  ditlerentiat  equation  lor  the  loaAtg  index  y  This  is  a  major  difference  betweeri  the 

present  theory  and  the  classical  theory  o(  plasticity.  The  loaning  irydex  now  should  be  treated  as  an 
Independecd  variable  governed  by  equation  (9).  {This  could  be  treated  numerically  and  irworporated 
in  usual  fime  element  programs  wfMre  equatiorts  (9)  can  be  solved  independently  at  each  integration 
increment. 


The  One-Oimeosional  Problem  lanoth  and  Time  Scales 

In  order  to  lustrate  the  impiicattons  of  the  theory  outlined  in  the  previous  section  to  ft#  shear  banding 
problem,  vie  cortsider  the  simple  shearing  of  an  infinite  blodt  in  the  x  dirediott  We  assume  that  the 
state  ol  stress  is  pure  shear  such  that  S,  -  -  t(0  are  the  only  nonzero  compooenls  of  8.  Similarty, 

ty,mf2  are  the  only  nonzero  components  o(  0.  Moreover,  yve  assume  isotropic  hardening  with 
no  back  sbass.  Then  it  can  be  shoem  that  equations  (1)-(3)  reduce  to 

il.pv  ,  t.<5(Y-Tn,  <”) 

3y 

where  G  s  the  elastic  shear  moddus,  /  is  the  plastic  strain  rate,  and  v  Is  the  velodly  field  in  Ihe  x 
direction  Various  classes  of  material  behavior  are  considered  below. 


Eiastic-Ptastc  Materials 


a)  No  gradtonts  C,  >  C,  ■  0 

For  this  case  we  drop  the  dependence  on  y  from  equation  (5)  and  obtain  •  i/H.  where  H  ■  9  k/^/. 
Tr.in  equaions  (11)  (with  y  -  v,)  yield  a  rather  familiar  partial  differential  equation 


-t  OH 

■ 


(12) 


where  the  subscript  y  indicates  partial  dertvaUve.  The  type  of  this  dHferentlal  equaiton  is  governed  by 


me  sign  ol  H.  If  H  >  0  the  equation  Is  hypertjoic  and  traveding  waves  with  a  speed  equal  to  c  exist. 
However,  when  H  <  0  the  equation  becomes  elliptic  and  wave  trappirtg  occurs  {19J.  In  this  respect  the 
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problem  Is  'Hi-posed'  due  to  the  charoe  ol  type  occurring  in  pte  governing  differential  equation  as  the 
material  enter  the  sotienirrg  regime  (H  <  0). 

b)  With  gradients  C,  >  oonstard.  C,  >  oorrstant. 

The  problem  of  'Hl-posedness'  is  eliminated  by  introdudng  higher  order  gradients.  This  can  be  seen 
by  considering  equation  (5)  (with  no  dependence  on  ^  which  upon  neglectlrrg  elasticity  from  (it), 
(y-yO  and  combining  with  (ii),  gnes 

pu^Hu„-C,u^-2C,u^„  (13) 

with  u  denoting  the  dispiaoement.  Now  the  character  ol  this  iWerentlal  equaHon  is  Irxlependent  of  H 
and  is  determined  by  C,  which  is  a  constant.  Therefore,  the  system  does  not  change  type  Moreover, 
equation  (13)  can  be  nondimensionaiized  by  introducing  the  lolowlng  ‘natural*  time  (t\)  and  lertgth  {/) 
scales 


;  p^i  .  ;  f'l 

H*  ■  IHl  ' 

Then  with  the  nondimensional  variables 

(14) 

f*-- 

II  n 

(15) 

equation  (13)  becomes 

■  ^rrr-r  ' 

(16) 

This  suggests  that  the  width  of  the  shear  bar>d  (w)  is  proportional  to  /,  t.e.. 


Wtien  H  >  0.  it  is  shown  in  [3]  that  the  homogeneous  state  is  stable  and.  therefore,  shear  bands  do  not 
initiate.  They  do,  however,  when  H  S  O',  i  e  according  to  (17)  at  H  •  0,  w  and  as  (H|  (in 
the  softening  regime)  w  -<  0.  corresponding  fo  severe  localization  This  is  consistent  with  ihe 
quantitative  result  given  in  (I2.t3j  where  equation  (16)  is  solved  numerically  with  the  inerba  term  on 
the  left  hand  side  set  equal  to  zero.  The  result  tor  70/30  brass  is  shown  in  Figure  1  and  compared  with 
the  experimental  ones  obtained  recently  by  Joshi  et  al  [20]  Note  that  such  nondimensionalizatlon 
canrrot  be  performed  on  equation  (12)  tor  which  a  "naturaJ*  letrgth  scale  cannot  be  defined 

£(asfic-viscop(astic  Materials 

Tor  this  case  we  drop  the  gradient  terms  from  equation  (5)  and  assume  that  (5),  can  be  inverted  to 
lormally  read 

y'’f(t,y^.  (’8) 


Then  equation  (18)  can  be  substrtuted  into  (it),  which  combined  with  (It),  yields 
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(19) 

or  p 

This  partial  differential  equation  Is  always  hyperboUc  sirtoe  the  speed  ot  the  propagating  wave  c  Is  real 
and  equal  to  the  ^>eed  ot  the  elastic  wave  (regardtess  of  tie  sign  of  the  hardening  modulus  This 
result  tor  elasttoviscoplastic  materials  is  dfffersnt  than  that  for  the  elastic-plastic  materials  given  by  (12) 
where  the  wave  speed  arxl  the  character  o<  the  govemirtg  ditferential  equation  Is  determined  by  the  sign 
of  H.  Moreover,  h  this  case  one  can  also  define  ’natural*  length  and  time  scales  as 

which  with  the  aid  of  the  nondimensional  variables  (15)  can  reduce  equation  (t  9)  Into  the  form 


where  vw/c.  This  nondimensionalization  implies  that  the  band  width  Is  proportion^  to  lor  In  view  of 
(20)  and  (t9)„ 

- - ^==-  (22) 

[atdrUtSp 

In  the  static  case,  however,  equation  (it),  implies  that  t  is  homogeneous  In  y  and,  therefore,  one  has 
only  equation  (18)  which  has  no  length  scale.  Hence,  higher  order  gradients  should  be  considered  and 
the  original  equadon  (5)  yields 

t(D*«<Y.f)-C,T^-<^y*  (23) 

whore  t(t)  Is  a  function  of  time  only  since  9t/3y  »  0  for  the  static  case.  Equation  (23)  has  been  solved 
In  [13]  numerically  where  it  is  shown  that  the  band  width  is  independent  of  the  mesh  size  and  depends 
only  on  the  values  of  C,  and  C,  as  can  be  seen  from  Figure  2. 


Effect  of  Plastic  I 


Numerical  Analysis 


For  the  numerical  analysis  considered  here,  we  neglect  elasticity  and  assume  that  the  material  is 
viscopiastlc  exhibiting  Wnematic/tsotropic  strain  harder^ng  and  softening  aooording  to  the  constitutive 
equations  (4),  (5)  and  (7).  This,  in  turn,  results  into  a  convenient  viscous  flow  FEM  formulation  (2t]  with 
a  nonlinear  secant  stiffness  nrtatrtx,  arid  eliminates  the  severe  numerical  limftations  associated  with  usual 
elasto-plastic  fonnolations.  in  fact,  equation  (7)  and  the  momentum  equation  (1),  with  the  inertia  term 
neglected,  along  with  the  usual  finite  element  formufatton  yield  [21) 

(K(X,V,TiIV.F,-  F..  (24) 


where  K  is  the  secant  stiffness  matrix  which  is  a  nonlinear  function  of  the  nodal  velocity  vector  V,  nodal 
position  vector  X  and  strain  gradients  r,.  F,  is  the  boundary  force  vector  and  F,  Is  the  body  totoe  vector 
arising  trom  the  back  stress  «.  The  nonlinear  dependence  on  V  arises  from  the  fact  that  the  viscosity 


p  is  nonlinearly  related  to  y  . 
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Equation  (24)  Is  sotvad  for  V  using  the  secant  method.  An  expHdt  Intogratioo  scheme  with  a  constant 
ttne  step  using  the  Newton's  forward  method  Is  used  to  integrate  the  velocity  and  the  evolution  equation 
for  the  back  stress  (4),.  This  integration  algorithm  along  with  the  secani  method  lo  handle  nortllneartties 
In  V  and  by  continuously  updating  the  geometry  at  each  time  inoament  to  handle  the  geometric 
nonllnearlties.  proved  to  be  very  efficient  and  produced  stable  results. 

For  the  numericai  analysis  we  use  the  following  power  taw  expression  to  model  the  isotropic  hardening 
and  softening 


«(t.y  )•», 


(25) 


where  t,  is  the  inifal  yield  stress.  m>Oandrt>Oarelhe  strain  rate  sensitivtiy  and  strain  hardening 
exponents  respectively,  v  <  0  is  the  'strain  softening'  exponent,  and  y»To  Ti  are  material 
oortstants  pt).  For  a  typical  stnrctural  steel  we  use  the  values  to=400l«^  MPa  ft  -  40''  MPA,  C  -  4.0, 


rt- 0.1.  m- 0.01,  V- -2.0,  To'2«t0  ’.  V  0  a«c''.  Ti'2.0  It  is  noted  that  lor  the  present  case 
the  gradient  effects  are  dropped  by  setting  C,  >  C, »  0. 


The  numerical  method  is  used  to  analyze  the  development  of  shear  bands  in  a  piane-strain  tensile  test. 
Due  to  symmetry,  only  one-quarter  ot  the  tensile  specimen  is  examined.  The  one-quarter  is  meshed 
into  350  rectanguiaf  elements  each  consisting  of  lour  cross-lrianguJar  efemerts.  The  specimen  is 
deformed  iri  tension  with  a  constant  cross-head  speed  of  V  ■  0.01  cnvSec.  A  length  to  width  aspect 
ratio  ot  3  is  used.  A  brief  summary  of  the  results  are  shown  in  Figures  3-7  which  are  sufficient  to 
Illustrate  the  effects  of  viscosity  and  plastic  spin  on  the  shear  bartds  development.  Figure  3  shows  the 
mesh  of  the  specimen  with  two  'cross'  shear  band  dearly  formed  A  doser  look  aA  the  bands  is  given 
in  Figure  4  where  only  one-quarter  of  the  specimen  is  shown,  it  can  be  seert  Irom  the  figure  that  the 
width  of  the  band  spans  over  at  least  two  elements.  The  correspondsig  strain  ooniours  are  shown  m 
Figure  5  iivdicating  V>e  severity  of  the  strain  m  the  band  and  the  tact  tuft  the  strain  IMd  varies  only 
across  the  band  wtih  constant  contours  along  It  The  eftecl  ot  the  plastic  spin  is  examined  by  varying 
the  value  of  (  [When  (-0  the  oorotatlana)  rate  given  by  (4),  reduces  to  the  usual  Jaumman  rate]. 
Figure  6  shows  P/P,  versus  U/L  for  various  values  of  a  where  P  is  load.  P,  is  yield  load,  U  is  elongation 
and  L  is  specimen  length.  The  sharp  drop  in  the  load-dispiacement  curve  for  U/L  >  0  2  corresponds 
to  the  development  of  the  shear  band.  The  band  seems  to  initiale  around  U/L  •  0.2.  It  seems  that  the 
parameter  a  has  no  influence  on  the  curve  until  the  band  becomes  severe  around  U/L  *  0  25.  This  is 
expeded  since  the  plastic  spin  influences  the  shearing  mode  which  (lor  the  present  problem  of  Initialty 
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stretching  model)  tiecomes  significant  only  v«hen  the  band  forms.  The  resui.  thus,  suggests  that  the 
plastic  spin  does  not  Influence  the  Initiation  of  shear  bands  from  a  stretching  mode  (see  also  Tvetgaard 
and  Qiessen  (22]),  It  may  homever  from  a  shearing  mode.  An  increase  In  the  degree  of  the  plastic  spin 
seems  to  decrease  'ductil%*  In  the  sense  that  the  value  of  the  average  strain  to  fln^  load  'collapse* 
decreases  as  a  ktaeases. 

Finalty,  Figure  7  shows  the  variation  of  the  kinetic  energy  E/E^  where  E,  ■  mV*/2.  m  Is  the  total  mass, 
and  V  Is  the  imposed  bounde ,  velocity .  it  can  be  seen  from  the  figure  that  the  Mnetic  energy  remains 
almost  constant  in  the  early  stages  of  deformation  hen  increases  steadily  urm  the  shear  band  initiales 
around  U/L  a  0.2,  resultirig  to  a  sharp  Increase  In  the  kinetic  energy.  Thereafter,  the  kinetic  energy 
reaches  another  'steady'  stte  once  the  barxf  Is  completefy  fomned.  Note  that  E/E,  1$  sBghtiy  less  than 
1 .0.  This  Is  expected  since  he  material  outside  the  band  moves  ifgidly  with  velocity  V. 
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FIgurr  1.  Shear  baoeKvidfb  versus  Strain 

in  7Q/30  Brass, 

A)  rolling  angle  of  45*, 

B)  rolling  ingle  up  0”. 


Figure  2.  EITecI  of  tnesh  size 

A)  Ay  •  0.1,  B)  Ay  .  05. 


Figure  3. 


Deformed  Mesh. 


Figure  4.  Deformed  meshes.  Figure  5.  Strain  Contours. 


Figure  0.  Effect  of  pbstic  spin  on 

Load-esieBion  curve. 


Figure  7.  Effea  of  plastic  spin  on 
Kinetic  energy. 
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